Crossover of the weighted mean fragment mass scaling in 2D brittle fragmentation 
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We performed vertical and horizontal sandwich 2D brittle fragmentation experiments. The 
weighted mean fragment mass was scaled using the multiplicity fi. The scaling exponent crossed over 
at log /j, c ~ —1.4. In the small p(<C /i c ) regime, the binomial multiplicative (BM) model was suitable 
and the fragment mass distribution obeyed log-normal form. However, in the large /j,(^$> /fc) regime, 
in which a clear power-law cumulative fragment mass distribution was observed, it was impossible 
to describe the scaling exponent using the BM model. We also found that the scaling exponent of 
the cumulative fragment mass distribution depended on the manner of impact (loading conditions): 
it was 0.5 in the vertical sandwich experiment, and approximately 1.0 in the horizontal sandwich 
experiment. 

PACS numbers: 46.50.+a, 62.20.Mk, 64.60.Ak 
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The origin of the power-law distribution in brittle frag- 
mentation is one of the best-examined problems in sta- 
tistical physics 0, Q- K has been examined in many 




of fragmentation transition and low-impact energy frag- 
mentation have been discussed 0, 0] . Due to the suc- 
cess of scaling theory with critical phenomena, it is nat- 
ural to consider the universality of critical behavior for 
various phenomena. Kun and Herrmann discussed the 
possibility of percolation universality using a point im- 
pacted granular solid model ,9j. They also investigated 
the universality of shell fragmentation ^(j- Astrom et 
al. proposed another universality law for LJ liquid and 
elastic beam models ^3|. Dimensional analyses of the 
exponent of the power-law distribution have also been 
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Previously, we conducted 2D brittle fragmentation ex- 
periments in which we applied a flat impact to one side 
of the specimen ^t|. This consisted of a vertical sand- 
wich procedure using glass tubes. We showed that the 
critical scaling differed from that of percolation transi- 
tion, and proposed a binomial multiplicative (or biased 
cascade) model for critical fragmentation. The binomial 
multiplicative (BM) model is very similar to the turbu- 
lent multifractal p model [Tsj . This implies the similarity 
between brittle fragmentation and turbulence by means 
of multifractality. However, the BM model included a fit- 
ting parameter that was fixed at a = 2/3, although the 
origin of this value was not clear. When a more realistic 
case was considered, the model predictions did not fit the 
experimental results 0] . The model results also did not 
follow the power-law fragment mass distribution; rather, 
they obeyed a log-normal distribution due to the central 
limit theorem. 
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Low-impact energy fragmentation measured in exper- 
iments that involved dropping a ID glass rod yielded 
log-normal distributions in the relatively low-impact en- 
ergy regime Q. The log-normal form has also been ob- 
served in the 31? numerical results of viscoelastic crys- 
tal fragmentation The first discussion of the log- 
normal distribution for a fragmentation process is found 
in Kolmogoloff 20] . It is not clear how the fragment 
mass distribution approaches the power-law form from 
the log-normal distribution. Do the fragments obey any 
other distributions before they reach the power-law form? 
The relation between the universal scaling law, the log- 
normal model, and multifractality is one of the most fre- 
quently discussed top ics, even in the turbulent energy 
cascade problem [21j. Since the brittle fragmentation 
phenomenon is very simple, it is very useful to investi- 
gate the origin of and path to the power-law form. 

In order to study this problem, we performed low- 
impact energy fragmentation experiments. In addition 
to the glass tube results we reported previously 0, 0] , 
we also analyze the results for glass plate samples, which 
correspond to a horizontal sandwich procedure. 

The experimental apparatus was very simple. Sam- 
ples were sandwiched between a stainless steel plate and 
a stainless steel stage. Then, a heavy brass weight was 
dropped along guide pole s. This experimental system 
was described in Ref. After fragmentation, all the 

fragments were collected and their masses were measured 
using an electronic balance. We broke 25 new glass 
plates. Fifteen were 30 mm x 30 mm x 0.1 mm in 
size, and ten were 60 mm x 60 mm x 0.1 mm in size. 
We set the measurement limit for the minimum mass at 
0.001 g, but only analyzed the data for fragments down 
to m m i n = 0.01 g. This m m i n value is same as that used 
in the glass tube experiments. The glass tubes and plates 
corresponded to vertical and horizontal sandwich proce- 
dures, respectively. 

Let us introduce a critical divergence of the weighted 
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FIG. 1: The log of the dimensionless weighted mean fragment 
mass [M2 1 (Miirimm)] as a function of the log of the pseudo- 
control-parameter fi. The scaling crossed over from a log- 
normal to a power-law distribution regime around log /i c ~ 
— 1.4. While a satisfied a < 1 in the small /i(<C /x c ) regime, it 
exceeded 1 in the /1 >/i c regime. This implies that the BM 
model is unsuitable for the large /1 regime. 



mean fragment mass, 



M 2 



Mim min 
where Mk and /x are written as 



(1) 



M k = 22 m k n(m), ^ = m min — . (2) 
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where m and n{m) denote the fragment mass and frag- 
ment number of the mass m, respectively. Note that 
the summation in Eq. (J2J) includes the largest fragment 
mass. The left-hand side of Eq. Q also includes the fac- 
tor wiZjjj, which was not considered in the previous defi- 
nition of a (Eq. (4) in Ref. 0]). This factor is a normal- 
ization term for the weighted mean fragment mass and 
gives a dimensionless value. It does not affect the value of 
the scaling exponent. The multiplicity parameter /x was 
first introduced by Campi as a pseudo-control-parameter 
to analyze nuclear fragmentation |22j |. It indicates the 
dimensionless normalized fragment number. 

The entire plot of log[M 2 /(-Mim m i n )] vs. log /i is shown 
in Fig. n] The figure shows that the scaling crosses over 
around (log/x c ~ —1.4). There are also two divergent 
points in Fig. ^ which are likely due to experimental 
failure, such as an oblique impact. However, we did not 
remove these points, since we do not have clear criteria 
to distinguish between success and failure. In the regime 
/i ^ /i c , the scaling exponent a can be described using 
the previously obtained value a — 0.84 < 1 17] . The 
higher-order weighted mean fragment mass exhibited a 
multi-scaling nature and its exponent agreed with the 
one predicted by the BM model. Therefore, we expect 



FIG. 2: Log-normal form of the fragment mass distribution 
in the low-impact energy regime (log/x = —2.55 data). The 
sample was 150-mm long glass tube. The inset shows a semi- 
log plot of the same distribution. The dashed line indicates 
an exponential-like tail. 



the fragment mass distribution to obey the log-normal 
form in this regime. Figure shows an integrated log- 
normal form of the cumulative fragment mass distribu- 
tion N(m) = J n{m!)dm! for a typical low-impact en- 
ergy fragmentation (150-mm long glass tube data with 
log/i = —2.55). The integrated log-normal function can 
be written as 



N(m) = A 



exp [-{log(m'/m)} 2 /2<7 



In] 



m 'V 27T<J fn 



dm' (3) 



where A, m, and cr;„ are parameters, which were taken 
as 0.24, 10.0, and 2.0 for the solid line in Fig. El respec- 
tively. We used = 20 as the cutoff scale. Since good 
agreement was obtained, the fragment mass distribution 
in fi <C fi c followed a log-normal distribution. 

While most of the data exhibited a log-normal form, 
there were a small number of fragments in the low-impact 
energy regime in general (e.g., raw curves in Fig.Efa)) so 
that it was difficult to establish the form of the distribu- 
tion directly. Therefore, we measured the weighted mean 
fragment mass using the moment Mk of the distribution 
to obtain sufficient evidence. However, we encountered 
problems when calculating the multiscaling exponent <7fc 
(defined as Mk+i / (Mkm min ) ~ H~ (Tk ) for the glass plate 
data due to the large fluctuations in M fe+1 /(M fe m m i n ). 
We did not obtain reliable estimates of tJfc for the glass 
plates, particularly in the large k region. Therefore, we 
focused only on M2/(Mim m i n ) scaling here. The scal- 
ing in the large k regime is obviously determined mainly 
by the largest fragment. This means that we used mean 
mass statistics instead of the largest mass statistics, in 
thia paper. 

For the glass tubes, some fragmentation results showed 
a power-law distribution in the relatively large (i regime 
|17j . In such regime, fi was close to /i c , i.e., the crossover 
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might have already occurred (Fig. 4(b) in Ref. 17]). Due 
to the dimensional restrictions of the experimental appa- 
ratus, we could only examine the small [i regime for glass 
tubes. Therefore, the clear crossover found in the glass 
plate fragmentation data has not been observed previ- 
ously. 

Another important characteristic is the power-law 
form of the cumulative fragment mass distribution for a 
fully fragmented state. It had different exponents for the 
tube and plate experiments. Figure shows the cumula- 
tive fragment mass distributions in the range m > 0.01 g 
for the glass plate samples. Figures|3Ja) and|HIb) give the 
low- and high-impact energy regime distributions, respec- 
tively. Each curve represents different imparted energy 
(dropping height of the weight) state. The cumulative 
distributions of well- fragmented events (Fig. Efb)) have 
a power-law portion N(m) ~ m - ' T_1 ) with an exponent 
t — 1 of about 1 . Some distributions in Fig. E^b) contain 
large fragments, so that the scaling regions are restricted 
to almost one order of magnitude; however, most por- 
tions of the distributions follow r — 1 = 1. The glass 
tube experiments had r — 1 = 0.5 This difference 

between the tubes and plates indicates that the exponent 
t depends on the fracturing method. The value r — 1 = 1 
does not concur with the value predicted by Hayakawa 
and Astrom et al, r - 1(= (d - l)/d) = 1/2 (for d = 2) 
0, . They considered the propagating and branching 
dynamics of the crack (or the failure wave). Therefore, 
in the horizontal sandwich fragmentation of glass plates, 
mechanisms other than crack dynamics might determine 
the value of r. Moreover, the boundary conditions are 
different between our horizontal sandwich experiments 
and the simulations of Astrom et al. By contrast, Be- 
hera et al. obtained a value of r — 1 ~ 1 in the highly 
fragmented state for a lateral impact disk fragmentation 
simulation |ll| |. This value agrees with our experiments, 
despite the difference in the loading conditions. Kadono 
discussed the energy balance and obtained the inequality 
1/2<t — 1 < 1 |6|. This inequality range is also close to 
our result. 

On the other hand, our distributions in the low-impact 
energy regime showed the remains of large fragments and 
were rather flat (Fig.[3fa)). This behavior resembles that 
of the integrated log-normal form, as described in Fig. 
□ Although all the curves in Figs. a) and fflb) corre- 
spond to different imparted energy states, we tried sum- 
ming up those. As a result, summed curves are shown 
in the insets; these more clearly indicate the integrated 
log-normal and power-law distributions. The solid curve 
in Fig. Ota) is same as that in Fig. [21 except for the cutoff 
scale moo = 3.5. 

The weighted mean fragment mass scaling of the glass 
plate samples only are shown in Fig.^c). Here, the tri- 
angles correspond to the low-impact energy regime dis- 
tributions (Fig. E£ a ))) an d the circles correspond to the 
high- impact energy regime distributions (Fig. Gib)). As 
expected, we confirmed a distinct separation between the 
two regimes using the weighted mean fragment mass and 
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FIG. 3: The cumulative fragment mass distribution of the 
glass plate samples in (a) the low-impact energy regime and 
(b) the high-energy-impact energy regime, (c) Scaling plot of 
log[M2/(.Mim m j n )] vs. log/i for all glass plate samples. The 
triangles correspond to the low-impact-energy cases (a), and 
the circles correspond to the high-impact-energy cases (b). 
The insets of (a) and (b) show the all summed curves. 
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the multiplicity. The value of a became large (a ~ 1.7) 
in the larger /i(^> /i c ) regime (Fig. Although the 
BM model can be applied to small a{< 1) values, it is 
inappropriate for large cr(> 1). Furthermore, other mod- 
els, such as the distributed and remaining cascade model, 
also break down for large a values [Tflj . The point [i c in- 
dicates the distribution crossover from the log-normal to 
the power-law. We cannot explain what happens in the 
large fj, regime at present. Perhaps the smallest limit of 
the splitting mass might appear above /i c , similar to an 
idea proposed by Matsushita and Sumida (23|. We as- 
sumed that the crossover point was universal, but there 
appears to be a slight difference between the points shown 
in Figs. H and 02c). More details and direct observa- 
tions of fragmentation are necessary to understand the 
crossover precisely. Theoretical studies are also required. 
In particular, an analysis in the vicinity of [i c would be 
interesting to see how the transition occurs. 

The crossover in Fig. ^ is reasonable from the view- 
point of the limit point. The limit point (log/x, 
log[M 2 /(M 1 m m i n )])=(0, 0) corresponds to the com- 
pletely fragmented state. In such a state, all fragments 
are the smallest unit size fragments. While it is extremely 
difficult to achieve such a state (i.e., fragment mass dis- 
tributions exhibit power-laws in general), it can exist as 
an ideal limit case. If the BM scaling stretches until 
log/i = in Fig.^ log[M 2 /(-Mi77i m i n )] never reaches the 
value 0. This is a nonphysical state. Therefore, it is 
natural that the crossover point corresponds to a certain 
value fi c . 

Astrom et al. recently proposed a generic fragment 
mass distribution form that was composed of a power- 
law portion and an exponential portion |l4j . The for- 
mer originates from the dynamics of crack branching and 
merging, and the latter results from the Poisson process. 
Their proposed form also applies to the low-impact en- 
ergy regime. The inset in Fig. |2 depicts a semi-log plot 
of the same N(m) distribution that was explained using 
a log- normal distribution. It shows a straight (i.e., expo- 
nential) tail, which suggests that the Astrom model may 
also be suitable. However, from the viewpoint of the mul- 
tiscaling nature of critical fragmentation, the BM model 
and log-normal distribution are more plausible. Diehl et 
al. obtained a similar coincidence between 2D explosive 
fragmentation simulations and the BM model ■ They 
also discussed the log-normal distribution form. 

Very recently, Wittel et al. reported the results of shell 
fragmentation experiments and simulations [Toj . They 



concluded that the impact fragmentation of shells showed 
a continuous transition, while the explosive one showed 
an abrupt transition. In our experiments, fragmentation 
seemed to occur suddenly. We could not obtain samples 
that only had visible macro-cracks, but did not split. 
A small amount of imparted energy cannot make brit- 
tle solids cleave. This might imply a "latent-heat-like 
behavior" . That is, the beginning of fragmentation re- 
quires a finite "latent energy" to generate macro-cracks. 
The splitting occurs abruptly and it proceeds according 
to the BM model statistics. We can observe critical scal- 
ing in the range [i > 0. However, we cannot discuss the 
scaling in the range /i < 0, since it corresponds to the un- 
fragmented state. The fragmentation transition of open 
2D objects involved in flat impacts is not yet understood 
very well in terms of the phase transition, and this is 
still an open question. Conversely, the transition from 
the log-normal to the power-law is characterized by the 
crossover of the weighted mean fragment mass scaling, as 
demonstrated above. 

Wittel et al. also revealed that the scaling exponent r 
is dependent on the loading conditions in numerical sim- 
ulations [ljj. While this was not consistent with their 
experimental results, it concurs with our findings qualita- 
tively if we consider the vertical and horizontal sandwich 
procedures to correspond to impact and explosive frag- 
mentation processes, respectively. Quantitatively, their 
values of r differed from ours slightly. This might result 
from the difference between an open 2D sample and a 
closed shell sample. 

In summary, we examined 2D brittle fragmentation us- 
ing experiments with glass tubes and glass plates. The 
exponent r had different values depending on the load- 
ing conditions, which consisted of either a horizontal or 
vertical sandwich impact to the 2D surface. Contrar- 
ily, the normalized weighted mean fragment mass scaling 
was universal and had a crossover point at which the 
fragment mass distribution changed from a log-normal 
to a power-law type. The results were consistent with 
other recent experiments and numerical simulations, but 
included new experimental findings about the relatively 
large /x weighted mean fragment mass scaling. 
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